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The great reduction in thermal conductivity of skutterudites upon filling the "void" sites with 
Rare Earth (RE) ions is key to their favorable thermoelectric properties but remains to be un- 
derstood. Using lattice dynamic models based on first principles calculations, we address the most 
popular microscopic mechanism, reduction via rattling ions. The model withstands inelastic neutron 
scattering and specific heat measurements, and refutes hypotheses of an anharmonic RE potential 
and of two distinct localized RE vibrations of disparate frequencies. It does indicate a strong hy- 
bridization between bare La vibrations and certain Sb-like phonon branches, suggesting anharmonic 
scattering by harmonic RE motions as an important mechanism for suppression of heat conductivity. 
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The discovery of new high performance thermoelectrics 
with complex crystal structures poses an important chal- 
lenge for theory, which, if met, could lead to the synthe- 
sis of materials with substantially higher thermoelectric 
performance. The key is to learn how the contradictory 
criteria — high electrical and low thermal conduction — 
are satisfied. One of the most promising groups of ma- 
terials are the filled skutterudites. Skutterudites (MPns: 
M=Co,Rh,Ir; Pn=P,As,Sb) can be described as a cubic 
lattice of M atoms with 3/4 of the cubes filled by nearly 
square Pn4 rings. They have reasonable electronic prop- 
erties for thermoelectric application, but ordinary values 
(«100 mW/cmK at r=300K) of thermal conductivity It 
was conjectured by Slack that a great reduction in ther- 
mal conductivity could be achieved by filling the empty 
cubes (voids) with rare gas atoms. The idea was that the 
fillers would "rattle" inside the voids, strongly scattering 
low frequency phonons but not charge carriers. It was 
subsequently found that filling with rare earths (La, Ce) 
reduces the thermal conductivity by a factor of more than 
five at room temperature. |l^;0| Essentially, regarding the 
thermal conductivity, these crystalline materials may be- 
have like glasses. The origin of this strange behavior, 
also seen in some other materials (clathrates, KBr-KCN), 
is basically not understood. The occasionally discussed 
possibility that such a large reduction of the lattice ther- 
mal conductivity (LTC) is due to scattering by electronic 
excitations can be rejected by various arguments based 
on experimental trends 0]. On the other hand, sorting 
out other mechanisms requires microscopic analysis. The 
most notable hypothesis relates to anomalous dynamics 
of RE atoms, namely strongly anharmonic rattling mo- 
tion. It has also been suggested, from experiment, that 
there exist two inherent RE vibrations (or perhaps two 
different two level systems (TLS's)) of significantly differ- 
ent frequencies implying strong anharmonicity of a dif- 
ferent nature than simple rattling. Here, we construct a 



microscopic model needed to address these issues, demon- 
strate that it describes the experiments, and use it to 
discuss possible scenarios for the LTC reduction. 

Our lattice dynamical model is based in part on first 
principles density functional calculations, performed by 
the linearized augmented planewave (LAPW) method 
Q , as described in Rcf. Q , including directly calculated 
forces generated by numerous (w 40) sets of atomic dis- 
placements small enough to represent the harmonic cou- 
pling constants along various directions. Previously, we 
B used the LDA to obtain Ag and A,j vibrational modes 
for CoSba within a frozen phonon direct method and 
combined these with results of infrared measurements 
to develop a valence force field lattice dynamical model. 
Here we start with the model for CoSba, and make a 
minimal number of additions to, and adjustments of pa- 
rameters needed to reproduce our new LDA results with 
a reasonable accuracy. However, we shall first address 
the shape of the RE potential well. 

Fig. |l| shows the LDA total energy as a function of RE 
displacements in (La,Ce)Fe4Sbi2: This was calculated for 
a perfect crystal corresponding to the skutterudite (bcc) 
lattice with displacements of the REs along a trigonal 
axis. Not only is a double well excluded, but the anhar- 
monicity is of an ordinarily weak nature. In fact, addition 
of a small quartic anharmonic term is sufficient to fit the 
results to well within O.lmRy over the entire (±0.5 A) 
range. The harmonic frequencies are 68 and 74 cm^^ for 
Ce and La in the corresponding materials and the anhar- 
monic shift is less than a few cm"-'^ at T up to lOOOK. 
(We refer to these frequencies as "bare" frequencies; we 
emphasize that they are neither crystalline normal mode 
frequencies nor the "single-ion" frequency of a RE with 
all other ions fixed including the other REs). The crys- 
tal symmetry imposes isotropy on the harmonic term but 
not on higher order terms. We computed the quartic pa- 
rameter in the La filled material with La displacement 



along a cubic axis to be 10.7mRy/bohr^ as compared 
with 6.27mRy/bohr'' for a trigonal axis. These are the 
extrema of the quartic tensor. 
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FIG. 1. Calculated LDA energy vs. RE displacement 
along a trigonal axis of -RiiFe4Sbi2. The least square fit 
parameters are a=28.82 (24.64) mRy/bohr^ and 6=6.27 
(5.50) mRy/bohr'' for La (Ce). 
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FIG. 2. La (bottom) and Sb (top) weighted vibrational 
density of states (WVDOS) for LaFe4Sbi2. The dashed 
line shows Sb WVDOS for CoSba which is similar to that 
for FeSba within the same dynamical model. The vertical 
line denotes the "bare" La frequency. The high frequency 
region where the Fe atoms dominate is not shown. 

Now we present the full dynamical model obtained 
from our LDA results, including the above. We added 
central force constants between the RE ion and its clos- 
est Fe and Sb atoms to the valence force field model ||] 
for CoSba. We found this simplest approximation to be 
quite sensible by comparison with LDA forces for various 
sets of La, Sb, and Fe displacements. Thus we directly 
get values of —0.8 and 1.64 (10^ dyn/cm) for La-Fe and 
La-Sb force constants, respectively. The RE bare fre- 
quency is then reproduced by construction. This inter- 
mediate (Z)-model, based on a CoSba model of ReL ra] 
and augmented by the above parameters to account for 
La filling, is useful in analyzing the direct lattice dynam- 
ics effect of the RE. But, there is also a considerable 
indirect effect, by changes in other force constants upon 



La filling and our final (jr)-model includes these, specifi- 
cally a 30% reduction of the two central intra-square force 
constants and a 10% reduction of the M-Sb central force 
constants. Below « 140 cm~^ both models give similar 
results. A test of the JF-model is to compare frequencies 
of Aig phonons with those calculated by the LDA for a 
Ce comound (137 and 157 cm~^ ): the ^-model gives 
141 and 162 cm"\ and the X-model 154 and 179 cm~^ 
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FIG. 3. Dispersions along F-H. The widths are the rel- 
ative La character. Note the break in scale at 80 cm~^ 

The phonon dispersions of the unfilled material have 
rather fiat Sb derived optic branches quite close to the 
bare frequency of 74 cm^^ of the RE ions as seen in 
figures 4 through 7 of reference B. For example, one 
observed IR (F„ symmetry) frequency in the unfilled 
structure is 78 cm^^. F„ symmetry allows for RE dis- 
placements. We find strong RE-Sh coupling that drives 
predominantly Sb modes up in frequency and predom- 
inantly RE modes down. Fig. shows that the Sb 
weighted VDOS is greatly affected by La-filling in the 
80-100 cm"-'^ region. With higher frequencies, or sym- 
metries incompatible with La motion, Sb-like peaks are 
also shifted up, as the La-Sb force constant is positive, 
but the shifts are smaller. This is seen in the atom- 
type (J) weighted vibrational density of states (Fig. 0), 
Gj{ll)) = J2ij \^i{j)\'^^i^ ~ f^i); where the sum on i is 
over wave vectors and phonon branches and j is over all 



type J atoms in the unit cell. Figs, ga and b show the 
character of modes in the filled material in the spectral 
region of interest. The dash lengths are proportional to 
|eLaP where e^a is the La-projection of the (normalized) 
polarization vector; e.g., e\^ = 1 would mean the mode 
involves only the La sublattice - this would then have 
to be at the "bare" frequency of 74 cm~^. Strong cou- 
pling of Sb and RE motion is evident for all wavelengths. 
A resonant interaction between a longitudinal acoustic 
branch and a low lying {RE dominated) optic branch is 
found. Fig. || shows the region of the acoustic spec- 
trum where resonant scattering by RE ions in disordered 
material can occur. IHl 
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FIG. 4. INS spectral difference between La and Ce filled 
compounds. Differences in the high frequency spectral re- 
gion where the Fe and Co atoms contribute heavily are 
not shown. Theoretical results are based on model T (see 
text) and experimental results are from Ref. bj 

Armed with this microscopic model, we can discuss sce- 
narios for LTC reduction and related experiments, partic- 
ularly inelastic neutron scattering (INS) and specific heat 
(for T < 45 K) for RECoFegSbia [|. A rather sharp peak 
at 50 cni^-'^ and a quite broad peak centered at about 100 
cm~^ were obtained in a difference spectrum, i.e., sub- 
tracting spectra for La and Ce filled samples and noting 
that the cross section is much larger for La than for Ce, 
to find RE vibrational contributions, n This led to a 
natural assertion of two distinct RE vibrations. 

Using LDA calculations and our model, we can check 
the main physical assumptions made in the earlier anal- 
ysis. H First, it was assumed that there is no difference 
in interatomic forces between La- and Ce-fillcd materials. 
By LDA calculation we find that the self force constants 
of the La and Ce do differ by 15%, but, calculating the 
neutron spectra from our model, we find that the La 
WVDOS is indeed rather close to the difference in the 
INS between the La- and Ce-fiUed materials: two ma- 
jor peaks, centered at 50 cm"-'^ and 100 cm~^, in the La 
weighted VDOS (Fig. ||), are clearly seen in the calcu- 



lated differential INS spectrum (Fig. ^ as well. These 
correspond well to the two peaks found by Keppens et 
al. g] (As usual, the experimental peaks are broader due 
to instrumental and intrinsic materials effects like disor- 
der and anharmonicity.) Most importantly, the calcu- 
lated two peak structure was obtained without an addi- 
tional localized La mode. Instead, the La spectral weight 
was transferred by hybridization from the main La peak 
(downshifted from the bare frequency of 78 cm"'^ to 50 
cm~^), to Sb modes at 100 cm~^. 

Ref. H presents more indications of two RE frequencies 
(and possibly a TLS), namely specific heat differences be- 
tween CoSba and RECoFe3Sbi2 and small unusual tem- 
perature dependent features of elastic constants. Clearly, 
a reconciliation of our calculations that do not produce 
two separate La phonons, and the specific heat experi- 
ments is needed. (We cannot address the temperature 
dependence of the elastic constants in our calculations; 
this alone, however, cannot be considered as compelling 
evidence for TLS's.) Here we present new, improved ex- 
perimental specific heat data that supersede the earlier 
data of Ref. |^. The two sets agree below r=20K, but 
considerably disagree at higher T . In terms of the previ- 
ous Einstein modeling of the data for T < 50 K, we find 
that two frequencies are still needed to give the specific 
heat difference between the unfilled and filled materials, 
but the strength of the higher frequency component is 
much less. The experiments were done on several com- 
pounds; CoSba, LaCoFe3Sbi2 and LaFe4Sbi2. The syn- 
thesis procedure was previously described O. Heat ca- 
pacity data from 2-300 K were taken with a commercial 
Quantum Design system. Results from this system are in 
good agreement with published values for standards like 
sapphire and copper. Some of the new results are shown 
in Fig. 1^ |1^ along with our calculations. The new ex- 
perimental heat capacity for the filled material is much 
closer to that of the unfilled material than the earlier 
experiment. 

The agreement with the present theoretical results for 
Gv is remarkable, especially considering that the model 
does not have the two additional localized high-frequency 
modes. Resolution of this seeming paradox lies in an 
implicit assumption, which lacking force constant infor- 
mation had to be made in Ref. B, i.e. that interatomic 
forces are transferable between the filled and unfilled ma- 
terials and that specific heat differences between the ma- 
terials are solely from i?£' vibrations. However, our LDA 
calculations show an important contribution to the Sb 
dynamics from La-Sb forces; besides, the intra-square Sb 
interactions themselves are decreased by approximately 
30% upon filling. While new rather localized La modes 
do appear upon filling in our analysis (at ~ 50 cm~^), 
and they are reflected in the increased specific heat at 
the low temperatures, there are two more changes which 
greatly affect the specific heat at low temperatures: first, 
as discussed, certain Sb modes are shifted to higher fre- 



quency, which decreases the specific heat; second, the 
softening of the Sb-Sb forces increases the specific heat. 
At still higher temperature (50-300 K) the lattice specific 
heat approaches the harmonic Dulong-Petit value: There 
is good agreement with the models, but those measure- 
ments are not shown here. 
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FIG. 5. Comparison of specific heat measurements with 
the harmonic force models for T < 50K. The models, 1 
and 2, for filled materials are X and T . The inset shows 
the large effect of La modes in the T<10K region: The 
lower two curves were calculated from the 50 to 70 cm~^ 
spectral region exclusively. 



So the origin of the LTC suppression is plainly more 
subtle than prior speculations about strongly anharmonic 
RE motions. Aspects of our results suggest scenarios 
for the role of the RE ions: (1) There is significant har- 
monic interaction between La and Sb in comparison to 
the inter-square Sb force constants and most likely this 
holds for the anharmonic interactions as well [|ll| . Thus 
the rare earth vibrations, which are not heat carriers due 
to their flat dispersion, do interact strongly with heat 
carrying phonons. Slack and Galginaitis |13] pointed out 
that the effect of Raman-like scattering of phonons by 
magnetic impurity levels in CdTe materials could reduce 
LTC and the reduction would be greater with increasing 
temperature. Perhaps a similar effect takes place here 
via the cubic anharmonic interaction between heat car- 
rying phonons and the 50 cm~^ predominantly La vibra- 
tions even for a perfect crystal. (2) Typically there is a 
significant lack of complete rare earth filling in samples. 
Harmonic resonant scattering of phonons by " impurity" 
La (or La vacancy) vibrations could give short relaxation 
times for frequencies near the La modes. These are in a 
region of high diffusivity, Gv^ in the models. Also, the re- 



duction in Sb-Sb interactions on filling implies force con- 
stant disorder in the Sb sublattice of partly filled samples 
(lattice specific heats of LaFe4Sbi2 and LaCoFe3Sbi2 are 
extremely similar in the T<50K region so the reduction 
in Sb-Sb interactions seems due to RE filling, not substi- 
tution of Co by Fe) . 

To summarize, we find that the REs are in a har- 
monic well up to large displacements, implying that the 
simplest rattling ion models are inapplicable to these 
filled skutterudites; we explain the two La peak spectrum 
by harmonic lattice dynamics as follows: Without the 
dynamical La-Sb interaction the pure La modes would be 
concentrated around 70 cm~^, with a substantial number 
of Sb modes at slightly higher frequencies. Strong La-Sb 
hybridization pushes these groups apart, to ~ 50 and ~ 
100 cm~^. We suggest, on the basis of our calculations, 
possible mechanisms for RE induced LTC reduction. 
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